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ABSTRACT 

We present the analysis of 3 hrs of a rapid time series of precise stellar radial velocity 
(RV) measurements (cr = 4.5 ms~^) of the cool Ap star /3 CrB. The integrated RV 
measurements spanning the wavelength interval 5000-6000 A show significant varia- 
tions (false alarm probability = 10^'"') with a period of 16.21 min (z/ = 1028.17 /iHz) 
and an amplitude of 3.54 ± 0.56 ms~^. The RV measured over a much narrower wave- 
length interval reveals one spectral feature at A6272.0 A pulsating with the same 16.21 
min period and an amplitude of 138 ± 23 ms~^. These observations establish /3 CrB 
to be a low-amplitude rapidly oscillating Ap star. 
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I INTRODUCTION 

Kurtz (1989) first suggested that the cool magnetic Ap star 
j3 CrB was a prime candidate to be a rapidly oscillating Ap 
(roAp) star since it had stellar properties similar to several 
known oscillating Ap stars (HR 1217, 33 Lib, and 7 Equ). 
The issue of whether P CrB is a pulsating star has important 
implications for the excitation mechanism in roAp stars as 
well as the possible existence of an instability strip for roAp 
stars. 

Several investigators have searched for pulsational vari- 
ations in P CrB using photometric measurements with null 
results (Heller & Kramer 1988; Kreidl 1991). Recently, 
Kochukhov et al. (2002, hereafter K02) presented radial ve- 
locity (RV) measurements for f3 CrB spanning the wave- 
length interval 6105 - 6190 A. Measurements of most lines 
were constant to a level of 20-30 ms~^; however, intriguing 
evidence for RV variations were found for one spectral fea- 
ture, Fe I A6165.4 A. This feature had an amplitude of 71 ± 

II ms"'^ and a period of 11.5 ± 0.5 min. Before we can place 
j3 CrB in the family of roAp stars the K02 result must be 
confirmed since the false alarm probability for the detection 
was rather high (= 0.016). Here we present our own precise 
radial velocity measurements for /3 CrB and show that it is 
indeed a rapidly oscillating Ap star, but not with the period 
found by K02. 



2 OBSERVATIONS 

Spectral observations of (3 CrB were made for 3 con- 
secutive hours on the night of 31 January 1998 (JD — 
2450844.903220) using the "2dcoude" echelle spectrograph 
of the Harlan J. Smith 2.7m telescope of McDonald Observa- 
tory (TuU et al. 1995). According to the ephemeris of Kurtz 
(1989) this corresponds to magnetic phase 0.11. One hun- 
dred fifty-five observations were made with exposure times 
of 50 sees which yielded a typical signal-to-noise ratio of 
150 per pixel. The description of the instrumental setup can 
be found in previous papers (Kanaan & Hatzes 1998; Mkr- 
tichian et al. 2003). Briefly, the wavelength coverage was 
4700 - 6700 A at a resolving power of 60,000. Precise stellar 
radial velocities were achieved via an iodine absorption cell 
placed in the optical light path. 



3 RESULTS 

3.1 Integrated RV Measurements 

To achieve the maximum RV measurement precision we first 
used the full spectral region covered by the iodine reference 
absorption lines to produce a "mean" RV value for all spec- 
tral lines. In doing so we excluded the two reddest spectral 
orders that covered the wavelength range 6100-6300 A. Here 
the iodine absorption lines are weak and this degrades the 
RV precision. To determine the relative RV for (3 CrB each 
spectral order was divided into 60 chunks so that we could 
model the spatial (and temporal) variations of the instru- 
mental profile (IP) which can introduce significant RV er- 
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Figure 1. (Top) The Lomb-Scargle (L-S) periodogram of the RV 
measurements of /3 CrB in the wavelcngtli interval 5000-6000 A. 
The L-S power (unitless) is a measure of statistical significance. 
(Bottom) The discrete Fourier transform (amplitude spectrum) 
of the same RV measurements. In this case amplitude is in units 
of ms~^ 

rors. In each chunk the RV shift was calculated along with 
the IP using the IP reconstruction procedure outlined in 
Valenti et al. (1996). The RV measurements from all chunks 
were then combined weighted by the inverse square of the 
RV standard deviation for each chunk. RV measurements of 
an individual chunk that differed by more than 4<t from the 
mean were rejected in the calculation of the mean. 

The RV measurements had a long-term trend of un- 
known origin. This could be due to rotational modulation 
caused by the abundance spots known to occur on Ap stars, 
long-period oscillations, or a stellar or sub-stellar compan- 
ion. Observations with a longer timebase are needed to an- 
swer this. Since we are primarily interested in searching 
for pulsational variations on time scales of roAp stars (5- 
15 min), a second order polynomial was fitted to the long- 
term variations and subtracted. Variations in the sky trans- 
parency during the course of the observations resulted in 
some spectra having significantly lower signal-to-noise ratio 
(S/N) than the mean {S/N w 150). Observations having 
S/N < 100 (fifteen data points) were eliminated from the 
time series analysis. 

Figure 1 shows both the Lomb-Scargle (L-S) peri- 
odogram (Lomb 1976; Scargle 1982) and the discrete Fourier 
transform (DFT) of the residual RV time series for /3 CrB. 
These show significant power at 1028.17 /iHz (P — 16.21 
min). When trying to detect a periodic signal in time se- 
ries data and assessing its statistical significance the Lomb- 
Scargle (L-S) periodogram is a more appropriate analysis 
tool over the DFT. In this case the L-S power is not a mear 



sure of the amplitude of the signal, but instead is related 
to the statistical significance of a peak. Many researchers in 
the field of oscillating stars are more accustomed to dealing 
with DFT. In this instance the power is directly related to 
the true amplitude of the periodic signal in the data and can 
also be used to assess the noise level of the data. The DFT, 
however, has limited use in quantifying the significance of a 
signal. We thus show both the L-S periodogram and DFT 
(lower panel of Figure 1). 

The L-S periodogram of the /3 CrB data has power more 
than 15 which is the first hint that we are dealing with a 
statistically significant signal. Using the Eq. 14 of Scargle 
(1982) this power results in a false alarm probability (FAP), 
i.e. the chance tiiat noise is creating the signal, of ~ 2x10^^. 

The power in a L-S periodogram only gives a rough es- 
timate of the false alarm probability of a signal. A more 
accurate determination is to use a bootstrap randomiza- 
tion technique (e.g. Kiirster et al. 1997). This technique 
randomly shuffles the RV measurements keeping the times 
fixed. A Lomb-Scargle periodogram is calculated for each 
"random" data set and these data re-shuffled. After a large 
number of these shuffles, the fraction of the periodograms 
having maximum power greater than the observed data peri- 
odogram over the frequency interval of interest (in this case 
f — 200-5000 /iHz) represents the false alarm probability 
(FAP), or the chance that random noise can produce the 
observed power in the periodogram. Using 2x10® shuffles 
we establish a FAP = 10^'' for the peak in Figure 1. 

Figure 2 shows the residual RV measurements after re- 
moval of the long-term trend phased to the 16.2 min period. 
The top panel is the raw data while the bottom panel rep- 
resents phase-binned {A(j> ^ 0.05) values. The amplitude of 
these RV variations is 3.54 ± 0.56 ms~^ (/f-amplitude) . 

3.2 Narrow-band Measurements 

As first shown by Kanaan & Hatzes (1998) the RV amplitude 
of individual spectral lines in roAp stars can differ by factors 
of 10-100. Later studies established that elements of Nd and 
Pr in some roAp stars had the highest pulsational RV am- 
plitudes (Savanov et al. 1999, Kochukhov & Ryabchikova 
2001; Kurtz et al. 2003). To search for spectral lines that 
that might have high amplitude RV variations the individ- 
ual wavelength chunks (each chunk typically contained at 
most 1-2 spectral lines) prior to combining these in a mean 
velocity were examined. Because high RV amplitude lines 
can also be found in the reddest orders (excluded from the 
integrated RV analysis) these were also used to search for 
RV variations). A periodogram analysis was performed on 
the RVs from the individual wavelength chunks. If a peri- 
odogram showed significant power (Lomb-Scargle power > 
7) a statistical assessment was made using the bootstrap 
randomization technique (50,000-200,000 shuffles). 

Only one chunk (A6271.7 - A6273.3A) had highly sig- 
nificant variations with the same period as found in the inte- 
grated RV data. The spectral region for this chunk is shown 
in Figure 3. The RV amplitude (K) for this chunk, 138 ± 
23 ms~^, was significantly higher than the value for the in- 
tegrated RV measurements. Figure 4 shows both the Lomb- 
Scargle periodogram and DFT of the RVs from this wave- 
length chunk. Using 200,000 shuffles we estabUshed FAP w 
1.5 X 10~® for this peak. The pulsational phase for this 
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Figure 2. (Top) The residual RV variations phased to a period 
of 16.2 min. (Bottom) The phased-binned averages {A(j> fa 0.05) 
for data shown in the top panel. In both panels data points are 
repeated for the second cycle. 



Figure 4. The Lomb-Scargle periodogram (top) and DFT (bot- 
tom) of the RV variations for the spectral region (A6271.7 - 
A6273.3 A) 



chunk is identical to within the errors of the pulsation phase 
calculated using the integrated RV measurements. 

Wc tentatively identify the spectral features in the high 
amplitude chunk as being a Cell A6272.026 A line that is 
blended in the blue wing by a weaker CrII A6271.870 A fea- 
ture. 
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Figure 3. Spectral region of the wavelength chunk (A6271.7 - 
A6273.3 A) showing high amplitude RV variations. 



4 DISCUSSION 

Our precise RV measurements for /3 CrB establish that this 
star is with high probability a rapidly oscillating Ap star 
with an amplitude of 3.54 ± 0.56 ms~^ and a period of 
16.2 ± 0.7 min. Although additional observations are needed 
to confirm our result, we believe that the signal we have 
detected is real for a number of reasons: 

• We derive a very low false alarm probability. 

The probability that this is a false signal due to noise is 

~ 10~'^ as determined through the bootstrap randomization 
procedure. This is a more rigorous way of determining the 
FAP than strictly using the Lomb-Scargle power. This result 
represents the combined mean velocity of all spectral lines 
spanning the wavelength interval 5000-6000 A. Any system- 
atic or instrumental error would have to affect all wavelength 
chunks used in the analysis in the same way which seems un- 
likely. 

• The variations are not due to a few outliers. 

We think it is unlikely that a few outliers are causing 
our signal. (We define an outlier as deviating significantly 
from the mean RV value and not from the best fit sine wave 
through the data.) If only a few outliers are driving the 



power in the periodogram then this would be evident in a 
high false alarm probability. 

For a real periodic signal in the presence of noise, the more 
data one accumulates the more significant the detection be- 
comes, even with the presence of a few occasional outliers. 
Thus the power in a Lomb-Scargle periodogram should in- 
crease with increasing number of data points. This is demon- 
strated in Figure 5 which shows the Lomb-Scargle power at 
v = 1028.17 /iHz (circles) as a function of the number of data 
points used in the periodogram. For comparison wc gener- 
ated a synthetic signal consisting of a sine wave with the 
same period and amplitude as found in our data and sam- 
pled in the same manner as the real data. Random noise was 
added with the same rms scatter as our measurements. The 
crosses in Figure 5 show that the L-S power of the fake data 
(with a signal present) increases with increasing number of 
data points in the same way. 

An iodine absorption cell was used to perform our RV 
measurements because it is designed to eliminate, or at least 
minimize instrumental effects. However, one could argue 
that there might be residual instrumental shifts not cor- 
rected by the reduction process, or that our instrumental 
profile modelling somehow introduces a false signal into the 
data. We do not believe this is the case. 

We investigated whether periodic instrumental shifts were 
present in our data using the wavelength solution calcu- 
lated from the iodine absorption lines. A Fourier analysis 
of these shifts showed the strongest period at 176 minutes, 
and smaller variations with periods of 20 and 24 minutes. 
All periods were significantly different from the one found 
in our RV data. Furthermore, we also performed a Fourier 
analysis of all IP parameters used in our RV determination. 
We found no variations at the period coinciding with the 
one found in the RV data. This demonstrates that using 
the iodine cell technique does an excellent job of excluding 
instrumental variations. 

• The same signal is found in the \6272 A feature. 

The spectral orders used to determine the integrated RV 
did not include the order covering the A6272 A feature. Thus 
this spectral line is not contributing to the signal found in 
our integrated RV measurements. However, a separate anal- 
ysis of this spectral chunk did find the same 16.21 min pe- 
riod and with the same statistical significance as with the 
integrated RV meEisurements. Although this is not an inde- 
pendent confirmation (we are using the same data set, but 
different parts of the spectrum), this result and the above 
arguments argue strongly that the 16.21 min oscillations in 
/? CrB are real. 

Our RV measurements fail to confirm the results of K02 
on the Fe I A6165 A feature. The wavelength chunk contain- 
ing this spectral line had an rms scatter of 130 ms~^ (due 
to the weak iodine absorption lines in this wavelength re- 
gion). However, we had a factor of 3 more measurements 
than K02. Monte Carlo simulations indicate that we would 
have detected an 11.5 min period with an amplitude of 70 
ms"'^ and a FAP of 10~*. An amplitude of 50 ms~^ would 
have been detected with a FAP = 0.01. Thus we would have 
detected any variations in the Fe I A6165 A line if they were 
present at the same amplitude as reported by K02. 

We believe that the result of K02 on j3 CrB is spurious. 
The FAP is much too high (FAP = 0.016) and our measure- 
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Figure 5. The Lomb-Scargle power sA, v = 1028.17 /iHz as a 
function of number of data points used in the periodogram using 
real (circles) and fake data (crosses). The fake data was generated 
using a sine wave with the same amplitude (3.56 ms"'^) and pe- 
riod (= 16.21 min) as found in our data and sampled in the same 
manner. Random noise with the same scatter (cr = 4.5 ms~^) as 
our data was also added to the fake signal. 



merits have failed to find the presence of an 11.5 min period 
in either the integrated RV measurements, or in the analysis 
of the narrow wavelength chunks. More troubling is that the 
K02 period of 11.5 ±0.5 min found in /3CrB is very close to 
the values of periods (fti 11.7 min) found in several spectral 
lines in the roAp star 10 Acjl using the same instrument. 

The most likely explanation for the RV signal found in 
P CrB by K02 is that it is due to short term uncorrected in- 
strumental variations in the Gecko spectrograph used for the 
measurements. Circumstantial evidence for this comes from 
Matthews & Scott (1995) who reported an 11.1 min period 
in the RV variations in 7 Equ. This period is different from 
known pulsation modes and has never been confirmed by 
subsequent RV measurements from other investigations. Al- 
though the authors used a superimposed mercury emission 
line from an arc lamp to eliminate instrumental shifts, the 
11.1 min period could still be a residual instrumental effect 
due to the stellar light a calibration source having a slightly 
different optical path (not the case when using the iodine 
absorption cell). 

The RV measurements of K02, on the other hand, were 
made without a simultaneous wavelength calibration (either 
I2 absorption cell or simultaneous Th- Ar calibration) , rather 
the wavelength calibration used a Th-Ar exposure taken be- 
fore and after the time series used for the RV measurements. 
This cannot correct for any instrumental variations on time 
scales significantly shorter than the time between the two 
calibration measurements. Because K02 have failed to ex- 



elude an instrumental origin with a period of ~ 11.5 min for 
the RV variations in /3 CrB their results on this star (and 
possibly 10 Aql) should be considered suspect. An investi- 
gation of the short-term instrumental shifts of the Gecko 
spectrograph would be useful. 

Our RV meEisurements for /3 CrB seem to establish that 
this star is an roAp star with the lowest RV amplitude (3.5 
ms~^) and one of the longest periods (16.2 min). (3 CrB is 
also unique among roAp stars in that it shows high ampli- 
tude pulsatiorial in only one spectral feature, the A6271.9 A. 
This showed significant variations (FAP = 1.5xl0~^) with 
the same 16.2 min period and an amplitude of 138 ± 23 
m s^^ . We tentatively identify this feature as a blend of Ce II 
and CrII lines. Why /3 CrB has such a low RV amplitude, 
almost no high amplitude spectral lines, and a longer period 
mode compared to other roAp stars with similar spectral 
properties should provide clues as to the origin of the roAp 
phenomenon. We are currently planning observations of this 
star over a full rotation period not only to confirm our de- 
tection, but to search for rotationally modulated amplitude 
variations of the pulsations. 
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